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ABSTRACT 


This thesis investigates the Physical Layer performance of single-input single- 
output (SISO) wireless communications systems, as well as multi antenna techniques 
such as multiple-input single-output (MISO) and multiple-input multiple-output (MIMO) 
systems, the last two utilizing the Alamouti-based space-time block coding (STBC) 
technique. All cases are based on the IEEE 802.16-2004 standard with OFDM using 
different values of coding rates. International Telecommunications Union (ITU) channel 
models are selected for the wireless channel in the simulation process. The particular 
setting we are interested in is the case where partial Channel State Information (CSI) is 
fed back to the transmitter for optimal control on the transmission rate. The performance 
results of the simulated SISO, MISO and MIMO systems are compared among 


themselves. 
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EXECUTIVE SUMMARY 


To achieve effective use of bandwidth in mobile wireless communications, 
broadband techniques like Orthogonal Frequency Division Multiplexing (OFDM) have 
been developed in the last decades. Features of OFDM which make it attractive are its 
capabilities of easily adapting to all sorts of adverse conditions caused by multipath, 
interference and low Signal to Noise Ratio. Because of its advantages, OFDM is widely 
used in several standards adopted worldwide. IEEE° 802.11, IEEE” 802.16 and 
HiperMAN are some examples. These standards are also known by their commercial 


designations: Wi-Fi, WiMAX, and HiperMAN, respectively. 


The performance of OFDM can be further improved by data rate adaptation. In 
such techniques, the signal-to-noise ratio (SNR) at each subcarrier is measured and the 
information is fed back to the transmitter. This is normally designated as channel state 
information (CSI) feedback. By this technique we can also improve robustness against 


narrowband interference, since in this case only some of the subcarriers are affected. 


Other techniques used currently to improve the performance in wireless 
communications systems are based on multiple antennas on the transmitter and/or on the 
receiver. These schemes increase the capacity of a wireless link leading to higher data 
rates. In particular by single-input multiple-output (SIMO) configurations we obtain 
receiver diversity while transmit diversity is achieved in the multiple-input single-output 
(MISO) case. They enhance performance by space time coding (with MISO) or they 
mitigate fading (by SIMO). Furthermore, the utilization of multiple antennas permits the 
focus of the transmitted energy in specific directions by beam-forming, accomplishing 
spatial multiplexing of multiple users. The overall effects of multiple-input multiple- 
output (MIMO) can be summarized in terms of reduction of the bit error rate (BER), 
increase of system capacity and a more efficient use of the transmitted power. 
Furthermore the combination of OFDM with MIMO techniques turns out to be very 


attractive in terms of improved performance. 
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The goal of this thesis was to investigate the performance of MIMO OFDM 
communication systems with CSI feedback, in particular using the IEEE® 802.16-2004 
standard. Several steps were performed on the path to obtain this objective. Three 
systems were used to obtain the results presented in this thesis: SISO, 2x1 MISO and 2x2 
MIMO. All systems have the CSI feedback feature. The SISO system served as the 
baseline comparison for the remaining systems as those had additional features to achieve 
higher performances. Both the MISO and MIMO systems were implemented using the 


transmit diversity technique known as Alamouti’s scheme. 


In the simulations performed during this thesis work the International 
Telecommunications Union (ITU) channel models were selected, since they are the most 
frequently used power-delay profiles in simulation environments for modeling purposes. 
These models have three types of users: Indoor, Pedestrian and Vehicular. For each one 
of these types of users the ITU specifies two profiles of multipath: Profile A, with shorter 
time spread, replicates rural macro-cellular surroundings, while profile B reproduces an 
urban macro-cellular environment. For scenarios of micro-cells with radius less than 


500m, profile A is also suggested. 


The simulations were performed using SIMULINK® R2008a from MathWorks™. 
The results obtained in this thesis are presented in terms of BER curves and it was 
performed in several phases. First, the systems were developed and their performance 
was measured for each of the ITU user-channels profiles. From the results obtained 
several SNR thresholds vectors were defined to enable the systems’ utilization of the 
partial CSI feedback feature. This was named auto rate control mode. Under this 
configuration all systems were tested for ITU user-channels profiles of interest, in 
particular the indoor and pedestrian categories. A last phase of tests was conducted when 
shadowing effects on the channels were simulated to observe all systems’ response and 


performance. 


As expected, the simulations show that an increase in the number of antennas 
used in a wireless communications system enhances its performance and capacity. In 
particular the 2x1 MISO system outperformed the SISO system and was outperformed by 


the 2x2 MIMO system. The improvements in terms of SNR vary among the user- 
XVi 


channels profiles used in the simulations. In general, the MISO system presented a 2dB 
to 3dB average performance improvement compared to SISO system. The MIMO 
system also showed a performance improvement of the same magnitude when compared 
to the MISO system. In some cases the MIMO system reached peak performance 
improvements of 18dB when compared to the SISO system. As a result, the MIMO 
system capacity is higher than that of the MISO system and the capacity of this is higher 
than the SISO system. 
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I. INTRODUCTION 


A. BACKGROUND 


Eagerness to achieve higher goals has led mankind to an unprecedented 
technological level. Globalization may have started in the fifteenth century with the 
Europeans expanding their empires throughout the globe. Today, however, today we live 
the true meaning of globalization, particularly on information sharing in real time. Behind 
this high-performance global information distribution system lies the work of scholars, 
engineers and technicians, all fairly unknown to the common user. This global 
development in wireless communications has pushed to the limits of available resources 
such as bandwidth and capacity. The need for accommodating growing demands requires 
higher technical sophistication in wireless communications due to the environment in 


which they are employed. However, mankind wants more and more. 


Today’s wireless communications goals are narrower bandwidth, lower power 
consumption, higher data rates, and error-free data links. The physical world in which we 
live imposes some bounds to all these variables. In both civilian and military wireless 
communication applications, reliable data links are highly desirable. The need to operate 
such systems reliably in adverse conditions such as in dense urban environments or in 
hostile jamming scenarios led to the development of a number of advanced techniques at 
each layer of a communication system. In this thesis, we will focus particularly on the 


physical layer. 


Currently, the two big families of transmitted broadband signals through out the 
world are code division multiple access (CDMA) and orthogonal frequency division 
multiplexing (OFDM). The latter, in particular, is widely used in several standards 
adopted worldwide. IEEE® 802.11, IEEE 802.16° and HiperMAN (high-performance 
metropolitan area network) by the European Telecommunications Standards Institute 
(ETSI) are by far the most popular examples [1]. Often, these standards are known by 
their commercial designations: Wi-Fi, WiMAX and HiperMAN, respectively. The main 


reason for selecting OFDM in all the referenced standards is its capability of obtaining 
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better performance in multipath channel environments [2]. These circumstances are 
permanently found mostly in metropolitan and sub-urban areas, where the demand for 
broadband wireless communications is extremely high due to large population density. 
The performance of OFDM can be significantly improved by subcarrier data rate 
adaptation on the basis of the signal-to-noise ratio (SNR) at each subcarrier. In turn, this 
would improve robustness against narrowband interference, affecting only part of the 
subcarriers. Furthermore, OFDM allows effective implementation of single-frequency 


networks, which is particularly attractive for broadcasting purposes [2]. 


Conversely, the present-day use of multiple antennas on the transmitter and/or on 
the receiver has become a viable technique to increase the capacity of a wireless link [3]. 
OFDM systems obtain frequency diversity by using multicarrier modulation [1]. Systems 
with multiple antennas, in contrast, can create a number of independent channels leading 
to spatial diversity [1]. The main advantage of spatial diversity over time and frequency 
diversity is that it does not require additional resources of time and bandwidth [1], [3]. In 
addition, using multiple antennas we can focus the transmitted energy in specific 
directions by beam-forming, so that we can obtain spatial multiplexing of multiple users 
[1]. Particular configurations of interest are receiving diversity for the case of single- 
input multiple-output (SIMO) and transmit diversity for the case of multiple-input single- 
output (MISO). These configurations and the basic single-input single-output (SISO) do 
not allow spatial multiplexing. They enhance performance by space time diversity (with 
MISO) or they attenuate fading (by SIMO). In short, the advantages of multiple antennas 
systems can be listed as reducing the bit error rate (BER), increasing the system capacity, 
while enlarging the area of coverage and reducing the transmitted power [1]. By 
combining OFDM and multiple-input multiple-output MIMO schemes the overall 


wireless communication system performance is greatly enhanced. 
B. OBJECTIVES 


The goal of this thesis was to investigate the performance of MIMO OFDM 
communication systems with channel state information (CSI) feedback, in particular 


using the TEEE® 802.16-2004 standard. Several steps were performed on the path to 
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obtain this objective. Initially, a supplied SIMULINK® R2008a SISO OFDM model, 
from MathWorks™, was benchmarked. On the basis of this a SIMULINK® R2008a 
MISO OFDM model, MathWorks™, was modified and its performance measured. 
Finally, the MISO OFDM model has been extended to MIMO OFDM, and its 
performance measured. In all these models, a feature of partial CSI is fed back to the 
transmitter for optimal control on the transmission rate. The performance results of the 
simulated SISO, MISO and MIMO OFDM systems are compared among themselves. 
International Telecommunications Union (ITU) channel models were selected for the 
wireless channel in the simulation process, since they are the most frequently used 


power-delay profiles in simulation environments for modeling purposes [1]. 
C. RELATED WORK 


SIMO, MISO and in particular MIMO systems have been an exciting, active area 
of research in the past fifteen years because of their capability of increasing wireless 
communication capacity. It all started in 1994 when Paulraj and Kailath patented the 
“Increasing capacity in wireless broadcast systems using distributed 
transmission/directional reception [3],[4].” Telatar showed in [5] the benefits of using 
MIMO systems on additive Gaussian channel with and without fading. In the following 
year, Foschini introduced a “Layered space-time architecture for wireless communication 
in fading environments when using multiple antennas” in [6]. Later, together with Gans 
they presented some results in [7], on the limits of wireless communication in a fading 
environment using MIMO systems. In 1998 Tarokh et al. presented the performance 
criterion and construction of space-time block codes (STBC) in [8], pointing to high data 
rate wireless communications systems with multiple transmitting antennas. In that same 
year, Alamouti proposed a transmit diversity technique for MISO systems that provided 
the same diversity order as the maximal-ratio receiver combining (MRRC) of a 1x2 
SIMO system [9]. He also showed the possibility of implementing such a scheme for a 
2x2 MIMO system and pointed to a generalization of 2xM MIMO systems, where M 
is the number of receiving antennas. This scheme does not require any CSI feedback to 


the transmitter. His proposal became very popular due to its simple implementation and 


the performance enhancements obtained. This scheme was eventually adopted by the 
IEEE 802.16 Broadband Wireless Access Working Group as part of this standard and it is 
an option of implementation to improve the IEEE 802.16-2004 based communication 
systems [10]. Nowadays, MIMO OFDM is still a very active research area. Some studies 
are pointing to more elaborate techniques where CSI feedback is considered to further 
improve broadband wireless communication links. Linear pre-coding using CSI 
feedback, and per-antenna power and rate feedback to reach MIMO OFDM theoretical 


capacity were recently proposed in [11] and [12], respectively. 
D. THESIS ORGANIZATION 


This thesis is organized into five chapters. The present chapter is the thesis 
introduction. Chapter II provides discussion related to OFDM, multipath channels, 
multiple input and output channels, and Alamouti’s transmit diversity scheme. Chapter 
Il offers a description on the several models used to obtain the results presented in this 
thesis, which are explained and interpreted in Chapter IV. Finally, Chapter V presents a 


summary of the conducted work, the conclusions and suggestions for further research. 


YH. MULTIPLE INPUT MULTIPLE OUTPUT ORTHOGONAL 
FREQUENCY DIVISION MULTIPLEXING 


A. INTRODUCTION 


In this chapter, we will introduce the theoretical basis behind OFDM and channel 
models. The concepts of path loss, shadowing, multipath, and MIMO will be introduced. 
In particular, the MIMO approach will be derived from the basic SISO configuration. 
Furthermore, a description of the space-time coding known as Alamouti’s scheme will be 
presented. The necessary relations to its implementation in the transmitter and receiver 


for the 2x1 MISO and 2x2 MIMO configurations will also be established. 


B. ORTHOGONAL FREQUENCY DIVISION MULTIPLEXING 
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Figure 1. FDM Transmitter, after [13]. 


OFDM is one of the methods employed for efficient transmission of data blocks. 
Its concept is a refinement of frequency division multiplexing (FDM), where each carrier 
frequency is modulated with a separate information symbol, carrying one or more bits. A 
simplified representation of FDM transmitter is shown in Figure 1. The complex 
baseband representation of a particular time varying signal corresponding to the nth 


subcarrier with frequency f, can be written as 
“ane, (2.1) 


Summing all the signals we obtain an expression for the output of the FDM transmitter as 


s(t) = Six, Qe", (2.2) 


n=l 
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Figure 2. FDM Spectrum, after [2]. 


As depicted in Figure 2, in standard FDM a guard band between channels exists 
to avoid inter-carrier interference (ICI). It becomes clear that such a solution demands a 


great amount of bandwidth. To solve this particular disadvantage of FDM, in the 1960s 
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the concept of OFDM was proposed [2]. The original proposals pointed to the 
transmission of parallel data streams using FDM in which the channels no longer had a 
guard band in-between, but actually had an overlap. Figure 3 presents this concept. It is 
visible in comparing these two plots that a bandwidth gain of almost 50% is achievable. 
In creating an overlapping FDM, technique the problem of crosstalk between channels 


(sub-carriers) arises. 
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Figure 3. Overlapping FDM Spectrum, after [2]. 


2. Orthogonal Frequency Division Multiplexing 


To eliminate the drawback of crosstalk in an overlapped FDM spectrum, the 


subcarriers must be orthogonal to each other in the sense of a standard inner product. 


A baseband complex representation of the basic circuit for the OFDM transmitter 
with N_ subcarriers is shown in Figure 4. The complex baseband representation of a 


M-QAM signal corresponding to the nth subcarrier with frequency f, can be written as 


[14] 


a G)aber" 0<157 (2.3) 


where each symbol b, takes values from an M-QAM alphabet in the time interval T , 


which represents a symbol period. Then, the normalized output signal of the OFDM 


transmitter becomes the superposition of all subcarriers as 


N-l N-l 
s(t)=— Ys, ()= be" Ost s7. 2.4) 
n=0 n=0 
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Figure 4. Oscillator Based OFDM Transmitter, after [1]. 
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If this signal is sampled periodically with sampling interval T.=T/N and the 
spacing between subcarrier frequencies is established as 1/T (i.e. f, =n/T), then the 
normalized output signal of the OFDM transmitter calculated at samples t = kT. can be 
expressed as 

{Aol - ya 
(KT, )=— Di be N (2.5) 


n=0 


which by definition is the Inverse Discrete Fourier Transform (IDFT). 


It is difficult to assemble an OFDM transmitter with several modulators, one for 
each subcarrier, each containing its own oscillator, where the spacing between 
frequencies is critical to obtain the desired orthogonality. By reaching the result 
expressed in equation (2.5), it is clear that instead of constructing a traditional multi- 
oscillator based transmitter, it is far simpler to build such a system using an IDFT chip, 
generate the overall OFDM signal in baseband and digital format, and finally convert it to 
analog and translate it to radio-frequency (RF) before transmitting it to a channel. This 


conceptual OFDM transmitter is depicted in Figure 5. 
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Figure 5. IDFT Based OFDM Transmitter, after [1]. 


In Figure 6 we present the general spectrum of a single subcarrier of an OFDM 


signal, considering the use of squared pulses in the modulation process. Its magnitude 


spectrum is then of the form sinc [7 ( bay ie )] . As depicted in Figure 7, the spectrum of 


each subcarrier will have zero crossings at frequencies f,.=n/T with integers 
n=1,2,...,N—1, and the peak frequency at n=0. By forcing the subcarrier frequencies 
to be multiples of the symbol period, they will be orthogonal in this interval. It is visible 
that, by establishing the previous relations, at the maximum of each subcarrier all other 
subcarriers are zero. Such a relation will allow a demodulation process free of 
interference from all other subcarriers present in the OFDM signal. This is called ICI 
avoidance [2]. The OFDM demodulation process uses the Discrete Fourier Transform 


(DFT). A simplified block diagram is presented in Figure 8. 
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Figure 6. Sub-Carrier OFDM Spectrum, after [2]. 
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Figure 7. OFDM Signal Spectrum with Ten Sub-Carriers, after [2]. 
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Figure 8. DFT Based OFDM Receiver, after [15]. 


C. CHANNELS 


1. Additive White Gaussian Noise Channel 


The simplest wireless communication channel model is the additive white 


Gaussian noise (AWGN) channel. As depicted in Figure 9, the output of the channel r(t) 


is simply the sum of the input signal s(t) and the noise source n(t). 
r(t)=s(t)+n/(t). (2.6) 


The noise source has its main origin on the electronics of the receiver, but also 
from space noise. The basic thermal noise model assumes a power spectral density that is 


flat for all frequencies [16] 
N, 0 
Because of its uniform power spectral density noise power is designated as white 
noise. Its autocorrelation function is the inverse Fourier transform of the noise power 
spectral density, which is a delta function scaled by N,/2 [16] 


R,(t)=FNG (f)} = Po(e). (2.8) 


The power spectral density and its autocorrelation function are shown in Figure 
10 and Figure 11, respectively. This model is described statistically as a Gaussian 


process, for this reason it takes the name of additive white Gaussian noise. 
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Figure 9. Additive White Gaussian Noise Channel, after [15]. 
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Figure 10. White Noise Power Spectral Density, from [16]. 
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Figure 11. White Noise Autocorrelation Function, from [16]. 
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2. Linear Time Varying Channel 


The linear time varying channel model is an evolved AWGN channel, which 
includes physical characteristics such as multipath propagation. Starting with equation 


(2.6), the channel output s(t) becomes the convolution of the input with the time varying 


channel impulse response, plus the AWGN [15] 


r(t) =5(r)*h(r,t)+n(t)=| 


—0O 


+00 


h(z,t)s(t-r)dr+n(t) (2.9) 


where h(z,t) is the channel response at time ¢ in reaction to an impulse applied at time 


t—t. The linear time-variant channel with AWGN model is presented in Figure 12. 
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Figure 12. Linear Time Varying Channel with AWGN Model, after [15]. 


3. Large-Scale and Small-Scale Fading 


Some authors [16], [17] divide the wireless channel effects in two types of fading 
effects: large-scale and small-scale fading. Specifically, the effects that cause the 
received power to vary are caused by long, medium and short distance phenomenons. 
These are path loss, shadowing and fading, respectively [1]. The first two are considered 


large-scale fading and the last small-scale fading. 
a. Path Loss 


The power obtained at a receiver antenna, placed at distance d from the 


transmitter, is obtained by application of Friis’ free space loss equation [17] 
14 


_ GGA’ 


2 





P.(d) P (2.10) 


: (42d ) 


where P is the transmitted power, G, is the transmitter antenna gain, G. is the receiver 


antenna gain, and 4 is the wavelength. From equation (2.10) the path loss in dB can 


defined as 





2 
SSt | (2.11) 


(42d) 


In order to accommodate a number of environments causing propagation 


PL(d)=10log (2) = “10 


Tr 


losses we consider the general model [1] 
d Y 
P(d)= RPL( a) (2.12) 


where P.(d) is the received power at the receiver placed at a distance d from the 
transmitter, PL(d,) is the measured path loss at a reference distance d,, and y is the 
path loss exponential. Values for this parameter are presented in Table 1 and are 
dependent on the propagation environment. From equation (2.12) the empirical path loss 
in dB can be expressed as 


Plgg(d) =10409{ 2.) =[P1(4)] +10709{ © (2.13) 


r 0 
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Table 1. Path Loss Exponent vs. Environment, from [17]. 

















Environment Path Loss Exponent, y 
Free space 2 
Urban area cellular radio 2.7 to 3.5 
Shadowed urban cellular radio 3 to5 
In building line-of-sight 1.6 to 1.8 
Obstructed building 4 to6 
Obstructed in factories 2 to3 


b. Shadowing 


A further factor to consider on the received power at the receiver is the 
loss caused by random obstructions. These are normally caused by buildings, trees or 
vehicles that can appear between a transmitter and a mobile receiver radio-wave 
propagation path. Such effect is called shadowing [1]. From equation (2.12) we can 


incorporate the shadowing effect as a random variable 
d Y 
P (d) RPL d,)x{ 2) (2.14) 


where y is a sample of the shadowing random process. It is clear that the power at the 


receiver is now modeled in a way that original path loss is an expected received power 
value (or mean), and the shadowing is a random alteration around the expected value. 


Normally, the shadowing y is modeled as a lognormal random variable as [1] 
7 =10"", where x~ N(0,07) (2.15) 


in which the shadowing standard deviation o, of the Gaussian distribution is expressed 
in dB. Typically this parameter is of the order of 6 —12dB [1]. 
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From equation (2.14) the empirical path loss plus shadowing in dB can be 


expressed as 


PL,,,, (4B) = rotog{ =| PL(d,) |+ 1 og{ +[X] (2.16) 


rt 0 


where | PL(d,) | is as defined before and |X] is a zero mean Gaussian random variable 


with standard deviation o, expressed in dB. 


c. Fading 


Fading or specifically small-scale fading is influenced by a number of 
factors such as multipath propagation, speed of transmitter and/or receiver, speed of 


surrounding objects, and bandwidth of transmitted signal [17]. 


The relative transmitter or receiver velocities will introduce an apparent 
frequency change on the transmitted signal. This effect is designated as Doppler shift (or 


spread) and is expressed as [17] 


= vf. cosé (2.17) 


d 
Cc 


v is the relative speed between transmitter and receiver, f. is the transmitted signal 


carrier frequency, @ is the angle between the direction of motion and propagation, and c 
is the speed of light. A parameter which can be defined from the Doppler shift is the 


coherence time and is expressed as [1] 


LAS, (2.18) 


From this relation it is understandable that if the Doppler shift is large the channel will 
change quicker than if the Doppler is small. This means that the Doppler plays an 


important role on the amount of time that we can consider the channel unaltered. 


The multipath propagation is characterized as the sum of several received 


signals due to reflection on the propagation path. Some multipath channels will include a 


17 


line of sight (LOS) component and, in others, this component is negligible. In the 
literature statistical channel models for both cases (with and without LOS) are provided. 
These models are designated as Ricean fading channel and Rayleigh fading channel, 
respectively. Figure 13 presents the multipath channel with LOS. From this figure, it is 


clear that the received signal r(t) will be composed of the original signal s(t), plus 


delayed copies (two in this case) of the original signal 


r(t)=a,(t)s(t—7,)+a,(t)s(t-7,)+a,(t)s(t—-7,). (2.19) 





s(t) Transmit Receive rt) 














Figure 13. Multipath Channel with LOS, after [18]. 


In general the received signal can be considered a sum of attenuated, time 
delayed, phase shifted replicas of the original transmitted signal. The baseband impulse 


response of a multipath channel can be expressed as [17] 


a [i(2zf.r;(t)+9;(z,1)) | 
h,(t,t) = > a,(z,t)e 6(r-7,(t)) (2.20) 
i=0 


where we define excess delay as the delay between the first arriving multipath component 


T,=0 and the ith multipath component. The propagation delay is neglected when 


plotting multipath channel impulse responses. a,(z,t) and z,(t) are the real amplitudes 
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and excess delays of the ith multipath component at time ¢, respectively. The term 
2x f.t,(t)+9,(7,t) represents the phase shift due to propagation and also reflections or 


scattering in the path. 


Based on the power spread P(r) we can define the mean and the root 


mean square (RMS) value of the time delay as [17] 


ar LP (ayn 








jp pee 2.21 

ee Pe) re 

Tas =? (7) (2.22) 
where 

= Yar Pee. 

Pt Sua ; 22d 

Le? DP) ie 

Now we can define the channel coherence bandwidth B. as [1] 

B » d re a (2.24) 


Cc 


Sige, OF 


max 


The channel coherence bandwidth B. gives an approximation on how 


wide is the distance between two frequencies that still have a correlated frequency 
response. That is, it defines the range of frequencies where a specific channel is regarded 


as having equal gain, thus the designation frequency flat. 
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Table 2. Broadband Fading Parameters vs. OFDM Design Impact, after [1]. 
Parameter If “Large”? If “Small”? Design Impact 

Ift>T, 

Delay spread, Ifrt<«T, The larger the z to the 7, the 
frequency 

a frequency flat more severe the ISI 
selective 

1 If i >T, Guideline to subcarrier width 

Coherence If —<T, B. 

bandwidth, B. B,. = B,/10 and number of 

B frequency 

a frequency flat : subcarriers L >10B,,/B. 
selective 
Doppler As f,/B,. becomes non- 
Sea If vf. >c,fast If of, <c, - 1B. 
uf negligible, subcarrier orthogonality 
i fading slow fading } ; 
c is compromised 

Coh T. small necessitates frequent 

sO eR ue If 7 >T, If T. <T, fast 

time, channel estimation and limits T 

T slow fading fading 


Cc 


but provides greater time diversity 


All the parameters presented in this section are important when designing 


an OFDM system, since the channel characteristics will have an impact on the system 


performance. In Table 2 a summary of these parameters versus its design impact is 


presented. In this table, T is the symbol time and L is the number of subcarriers in an 


OFDM signal. Figure 14 shows the impact of both large-scale and small-scale fading on 


the received signal versus the transmitter to receiver distance. 
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Figure 14. Combined Path Loss, Shadowing and Multipath Fading, from [19]. 


4. Single Input Single Output Channel 


The SISO channel model, depicted in Figure 15, is the classical communication 
link between a single transmitting and single receiving antenna. As a starting point, let us 
consider equation (2.9) that was the mathematical representation of a received signal that 
had travel through a time varying channel with AWGN. From that starting point, for 
simplification purposes, we neglect the AWGN and consider that the channel is casual 


and its impulse response is finite with duration 1,,,,,. Then a SISO received signal 


becomes 


r(t)=s(t)*h(z,t) = i, h(z,t)s(t-c)dr. (2.25) 
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Figure 15. SISO Channel Model, after [1]. 
ay 


5 Single Input Multiple Output Channel 


The SIMO channel model has a single transmitting antenna and several receiving 
antennas i € {1, 2,...M a . The SIMO channel model is shown in Figure 16. In this case, 
each receiving antenna receives a signal that is the convolution between the transmitted 


signal and the ith channel impulse response h, (z,t) 


r.(t)=s(t)*h,(z,t),i=1,2,...,Mp. (2.26) 


t 


If we define both the received signal and the impulse responses as vectors of 


dimensions M, x1 as follows 


r(t)=[7,(t),7,(t),--7%, (4) | (2.27) 

h(r,t)=[h,(r,t),y(v.t),--hy, (7.0) ] (2.28) 
then equation (2.26) becomes 

r(t)=h(z,t)*s(t). (2.29) 












































Figure 16. SIMO Channel Model, after [1]. 
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6. Multiple Input Single Output Channel 


The MISO channel is basically the opposite of the SIMO channel model. In this 


case we now have several transmitting antennas j € {le 2,....M ph and a single receiving 


antenna. Figure 17 represents such a model. In this case the received signal is the sum of 


each transmitted signal convolved with its respective path impulse response 


TOPOS OLLTCL! (2.30) 


Jal 


Let us define the transmitted signals and the impulse responses as vectors of 


dimensions M,.x1 and 1x M,,, respectively 


s(t)=[s,(1), 5,(t).-5y, ()] (2.31) 


h(z,t)=[h,(.t),M,(c.t),--Ay, (7.1) (2.32) 


then equation (2.30) can be expressed as 


r(t)=h(z,t) *s(r). (2.33) 





























ii 


Figure 17. MISO Channel Model, after [1]. 
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ep Multiple Input Multiple Output Channel 

In the MIMO channel model, as depicted in Figure 18, we have several 
transmitting antennas j € ae 2,....M - and several receiving antennas i € {1, 2,....M s In 
this case, the received signal r.(t) in the ith receiving antenna is the sum of all 
transmitted signals throughout all transmitting antennas M, with all possible impulse 


responses h, ,(7,1) 
r(t)=>05,(t)*h, ;(2.t),i=1,2,---,Mg. (2.34) 


If we define all the possible channel impulse responses h, , (z,t) , between any jth 


transmitting antenna and any ith receiving antenna, in a matrix of dimensions M, x M, 


as 


H(z,t)= (2.35) 








| tu, (7>4) hyo (T.t) ooo Puy, .a, (752) | 


and s(t) as defined in equation (2.31), then the output of the MIMO channel can be 


expressed as 
r(t)=H(z,r) *s(t) (2.36) 


where r(t) is as defined in equation (2.27). 
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Figure 18. MIMO Channel Model, after [20]. 























8. Simulation of MISO, SIMO and MIMO Channels with Multipath 


A multipath channel can be simulated as the sum of tapped delayed lines 
originated form a single input. Each line incorporates complex exponential time variant 
coefficients. This model is shown in Figure 19 and generically implements equation 


(2.20). 


A MIMO 2x2 channel model with multipath, as depicted in Figure 20, can be 


expressed as [3] 


e® 0 | 1 pe ll A(z.) 0 1 p,|le* 0 
H(s)=| 0 ae 1 l 0 ale 1 I 0 al (2.37) 


where h, (z,t) and h, (z,t) represent the multipath tapped delay line components and are 
written in the form 


N-1 


h,(t,t) =h,(0,t)6(7)+ D_h, (i,t) 0(¢-7,), fore =1,2. (2.38) 


c 
i=l 


The phases ¢, with i=1,2,3,4, are chosen randomly and account for antenna position 


and the signal angle of arrival or angle of departure. The coefficients , and p, account 


2 


for the level of correlation between the antennas at transmitter and receiver. The tapped 


delay line components h, (z,t) are independent complex Gaussian random variables with 


delay z, and phase shift @,. 
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Figure 19. Multipath Channel Model, after [17]. 


The MISO and SIMO channel models can be constructed from a MIMO channel 


model by leaving a single output and a single input connected, respectively [3]. 
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Figure 20. MIMO 2x2 Channel Model with Multipath, after [3]. 


D. ALAMOUTYS SCHEME 


L, Alamouti’s Scheme in MISO 2x1 Configuration 


Consider a MISO system with two transmitting antennas and one receiving 
antenna. Define the complex baseband signal from the first transmitting antenna (tx 


antenna 1) as s,, and the signal from the second transmitting antenna (tx antenna 2) as 
s,. In the following symbol period however, let antenna one transmit —s, and antenna 


two transmit s;, where * represents the complex conjugate of the symbol. The space- 


time encoding described is summarized in Table 3, where T is the symbol duration. The 


complete MISO 2x1 system using Alamouti’s scheme is shown in Figure 21. 


ZT 


Table 3. | Encoding and Transmission Sequence of Alamouti’s Scheme, after [9]. 
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Figure 21. MISO 2x1 System using Alamouti’s Scheme, after [9]. 






































Now let us consider that the channels between the two transmitting antennas and 
the receiving antenna have flat fading behavior during two consecutive symbols, and that 


they can be modeled as complex exponentials. Then the channels can be expressed as [9] 
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h,(t)=h,(¢+7T)=h, =ae" 


. (2.39) 
h,(t)=h,(¢+7T)=h, =a,e"" 
and the received signal at time f and t+T can be written as [9] 
r=1r,(t)=r(t)=hs,+h,s,+n 
1 He ) ( ) “1 2-2. 1 (2.40) 


tr, =1,(t)=r(t +7) =—hs; +h,s; +n, 


where n, and n, are complex additive random variables that account for noise at the 


receiver. 


The combiner receives the received signals and the channels estimates and 


constructs the following relations [9] 


s, =hyr, +h,r, (2.41) 
8, = hor, —hy, 
which will be provided to the maximum likelihood detector. If we substitute equations 
(2.39) and (2.40) into the combiner equations (2.41), we can further expand these 
relations as [9] 


S, =(a; +a), +hjn, +h,n; ei) 


_ 2: 2: * * . 
S$, = (a; is Js, —h,n, +h,n, 


Finally, the maximum likelihood detector will choose the symbol s,, measuring 


the Euclidian distance between the received symbol and the symbols used in the 


constellation of the transmitted signal, and choosing the closest [9] 


(a; +a; -1)Is, : +d’ (5,,s;) <(@; +a; -1)}s,| +d°(8,,s,), Wi ee (2.43) 





The outputs of the maximum likelihood detector s, and $, are the estimates of the 


transmitted symbols s, and s,, respectively. 
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2. Alamouti’s Scheme in MIMO 2x2 Configuration 


The MIMO 2x2 Alamouti’s scheme uses the same transmitting relations as the 


MISO 2x1 scheme. However, since we now have two antennas on the receiver, instead of 


one, the channel and the received signal will double when compared to the previous 


configuration. This system is depicted in Figure 22. The channels to antennas relations 


and received signals notation is presented in Table 4 and Table 5, respectively. 
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Figure 22. 


30 


MIMO 2x2 System using Alamouti’s Scheme, after [9]. 


Table 4. Channels to Antennas Relations of Alamouti’s MIMO 2x2 Scheme, after [9]. 





Rx Antenna 1 Rx Antenna 2 
Tx Antenna 1 h, =a,e"" h, =a,e” 
Tx Antenna 2 h, =a,e” h, =a,e" 


Table 5. Received Signals Notation of Alamouti’s MIMO 2x2 Scheme, after [9]. 





Rx Antenna 1 Rx Antenna 2 
Time ¢ r ie 
Time t+T f, i 


Now we can write the received signals presented in Table 5 as [9] 


r, =h,s, +h,s, +n, 
r, =—h,s, +h,s, +n, (2.44) 
r, =h,s, +h,s, +n, 


r, =—h,s, +h,s, +n, 


where n, with k =1,2,3,4, are complex additive random variables that account for noise 


at the receiver. The combiner takes the received signals and the channel estimates and 


assembles the following relations [9] 


S, =hyy, +h,r, + hjr, +h,r, (2.45) 
8, =hjr—hy, + hjr,—h,r, 


that can be further expanded, if we substitute in equations (2.44) and the relations 


expressed in Table 4, as [9] 
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= (a; +a, +a;+ a; )s, +hin, +h,n} +hin, +h,nj 

7 Bs ecine’ Fattn ois (2.46) 
8, = (a; 5 1 Oe Oy )s, —h,n, +hjn,—-h,n, +hjn, 
The maximum likelihood detector will then choose the symbol s, and s, using 


the following relations, respectively [9] 


2 





S. 


1 





2 2 2 2 
(a; +a,+a;,+a, -1) 


<(aj +a +a;+a;-1 





2 
Ss 


i 








2 2 2 2 
(a; Has te Pay -1) 


( 
) se (2.47) 
( 


The outputs of the maximum likelihood detector §, and §, are the estimates of the 


transmitted symbols s, and s,, respectively. 


E. SUMMARY 


In the present chapter, we have introduced OFDM modulation and channel 
characteristics important to OFDM systems design. Practical implementation of channel 
models was also introduced. Finally, Alamouti’s scheme for MISO 2x1 and MIMO 2x2 
systems was discussed. In the next chapter, we will describe the SISO, MISO and MIMO 
OFDM models used in simulations performed in this thesis to obtain the results presented 


in Chapter IV. Furthermore, a presentation on the ITU channel models will be provided. 
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HiIl. MODELS DESCRIPTION 


A. INTRODUCTION 


In the previous chapter, the basic theory on channel models and OFDM 
modulation was presented. In this chapter, a description of the models used to obtain the 
results later presented in this thesis will be provided. The basis of the work conducted 
during the thesis is a model provided by the MathWorks™ in the MATLAB® & 
SIMULINK® R2008a software package. Specifically, the IEEE® 802.16-2004 OFDM 
PHY Link, Including Space-Time Block Coding model was used. This base model can be 
found under Demos-Blocksets-Application Specific Examples. The base model was 
modified and improved on several blocks and several settings were changed to meet the 
research and simulations requirements. As a result of the work presented in this thesis 


three models were developed as described in the following sections. 
B. SISO OFDM MODEL 


The SISO OFDM model is shown in Figure 23 and its blocks are presented in the 


following sections. 
1. Forward Error Control and Modulator Bank 


This bank comprises seven FEC and modulator lines, one for each of the seven 
modulation schemes and the overall coding rate used in the IEEE® 802.16 standard: 
BPSK r=1/2; QPSK r=1/2; QPSK r=3/4; 16QAM r=1/2; 16QAM r=3/4; 
64QAM r=2/3; 64QAM r=3/4. 


a. BPSK r=1/2 FEC and Modulator 


The stream of random bits generated by the Data Source block is provided 
in blocks of bits designated unconstructed blocks. Several data source blocks are 
collected to form the data blocks provided to the modulator .The sizes of these blocks are 


presented in Table 6. The last column of this table presents the final value of the blocks 
33 


of data bits after appending the tail byte of zero. The blocks of bits are encoded by a 


convolution encoder of rate r=1/2 and generator polynomials g,=171,., and 
g, =133,,. No puncturing is executed. The data is then interleaved and finally 


modulated with BPSK modulation. 
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Figure 23. SISO OFDM Model. 
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Table 6. Data Block Size per Modulation Scheme. 














Modulation sehene Unconstructed —_ Block size before Block size after 
block size zero pad tail byte —_zero pad tail byte 
BPSK r=1/2 12 88 96 
QPSK r=1/2 24 184 192 
QPSK r=3/4 36 280 288 
16QAM r=1/2 48 376 384 
16QAM r=3/4 72 568 576 
64QAM r=2/3 96 760 768 
64QAM r=3/4 108 856 864 


b. Remaining FEC and Modulators 


In the remaining six FEC and modulators lines the stream of random bits 
generated by the Data Source is grouped in blocks and padded as described in the 
previous section. The blocks of data are afterwards applied to a Reed-Solomon encoder 
with the specifications presented in Table 7. Following that operation the data is encoded 
by a convolutional encoder with generator polynomials presented in the previous section. 
This encoder also provides puncturing in order to obtain the overall coding rates 
previously stated. In Table 8, the convolution encoder puncturing configuration is 
presented and its conceptual implementation is depicted in Figure 24. The data is then 
interleaved and finally modulated with one of the modulation schemes presented in the 
first column of Table 7. Finally, the last block yields the respective modulation scheme 
block, i.e., the QPSK’s lines have QPSK modulator blocks and the M-QAM lines have 
M-QAM modulator blocks. 
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Encoder 


[171,133] 





roo =1/2, P=[LI] 
loo =2/3, P=[L1,0,1] 

roo =3/4, P=[I,1,0,1,1,0] 

roc =5/6, P=[I,1,0,1,1,0,0,1,1,0] 





Figure 24. Convolution Encoder and Puncturing, after [18]. 


Table 7. Reed-Solomon Encoder Specifications, from [10]. 

















Modulation Overall Reed- Convolution 
scheme coding rate Solomoncode code rate 
QPSK r=1/2 1/2 (32,24) 2/3 
QPSK r=3/4 3/4 (40,36) 5/6 
16QAM r=1/2 1/2 (64,48) 2/3 
16QAM r=3/4 3/4 (80,78) 5/6 
64QAM r=2/3 2/3 (108,96) 3/4 
64QAM r=3/4 3/4 (120,108) 5/6 
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Table 8. Convolution Encoder Puncturing Configuration, from [10]. 














Code rates 
Rate 1/2 2/3 3/4 5/6 
fice 10 6 5 4 
X 1 10 101 10101 
Y 1 11 110 11010 
XY XY X VY, XVY,X,  X,VY,X.Y,X. 


2. OFDM Modulator 


In the IFFT Input Packing block the data from the modulators is first transformed 
from serial to parallel and then distributed through the sub-channels. At this point the 
data corresponding to the DC carrier and the Pilot Carriers is generated. The DC index is 
simply set to zero in order to decouple the data from the carrier. The pilots are all 
produced using the same pseudo-random sequence, BPSK modulated as produced by the 


Pilot Generator. The indexing for the IFFT input is shown in Figure 25. 


In the OFDM transmitter, first a short preamble is appended to the data provided 
by the IFFT Input Packing. The preamble, in all three models, is used for synchronization 
and to estimate the CSI. Then twenty-eight nulls in the lower indexes and twenty-seven 
nulls in the upper indexes are added to the original two-hundred-and-one indexes to 
provide frequency guard bands. The ordering of the data follows the modulo operation 


(256 carriers in this case) and it is shown in Figure 25. 


A gain block is used to compensate for the number of sub-carriers and normalize 


the transmitted power. Its expression is given by 





Gad, (| 56. > 18.1. (3.1) 


Li iced 
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Figure 25. Inverse Fast Fourier Transform Input Packing, from [18]. 


3. Multipath Fading Channel with Additive White Gaussian Noise 


This block can simulate three different channel configurations: AWGN channel 
only in which only noise is added to the signal; Frequency-flat fading with AWGN; 
Frequency-selective fading with AWGN. The simulations executed only the first and the 


last options, mainly. The fading mode and SNR in dB are defined by the user. 
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4. OFDM Receiver 


In the OFDM receiver the data received from the channel is reshaped from serial 
to parallel, following the removal of the cyclic prefix. The data is then transformed 
through a FFT operation and then rescaled as in the OFDM transmitter to compensate for 
the number of sub-carriers used. This gain is given by equation (3.2). The data is then 


extracted after the FFT following an ordering analogous to Figure 25. 


ie it 256 


: ‘FFT _ 


1 
G= =. 
V Lig Li sed 256 200 





= 0.07071. (3.2) 


oe Gain and Phase Compensator 


In this block the gain and phase of the received signal is corrected by comparing 
the original transmitted preamble with the received preamble. The calculated values of 
compensation are applied to the received signal after the DC component removal. The 


compensated signal is then sent to the data carriers’ extraction block. 
6. Data Carriers Extraction 


In the data carriers’ extraction block the pilot subcarriers are extracted and 
dumped. Recall that in this model the pilot subcarriers are not used for channel 
estimation, instead, the preamble is used for that effect. Finally, the data is converted 
from parallel to serial and sent to the demodulator and FEC bank and also to a 


constellation scope. 
7. Demodulator and Forward Error Control Bank 


As the FEC and modulator bank, this block is built with seven lines, one for each 


of the modulation schemes used. 
a. BPSK r=1/2 Demodulator and FEC 


In this line the data is demodulated, deinterleaved and decoded using a 


Viterbi decoder. Hard decision decoding is employed. The deinterleaver and Viterbi 
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decoder settings match those of the interleaver and convolution encoder used in the FEC 
and modulator bank, respectively. This block also modulates the demodulated data to be 
summed to the inverted received modulated data. This result is provided to the SNR 


estimation block for channel instantaneous SNR estimation. 
b. Remaining Demodulators and FEC 


The remaining six demodulators and FEC are implemented in the same 
form as the BPSK r=1/2 demodulator and FEC, with the addition of the punctured 
Reed-Solomon decoder. Again, the hard decision decoding is used in the Viterbi decoder. 
Also, the deinterleaver, Viterbi decoder and punctured Reed-Solomon decoder settings 
match those of the interleaver, convolution encoder and Reed-Solomon encoder used in 
the FEC and modulator bank, respectively. The settings used are those presented in Table 
7 and Table 8. For the remaining schemes, only the modulator and modulator blocks are 


replaced by other blocks with the respective modulation scheme. 
8. SNR Estimation 


In this block the instantaneous link SNR is calculated based on the estimated error 
between the transmitted and received symbols. This information will be used for adaptive 
rate control of the system. The root-mean-square of the column vectors of the signal is 
calculated and then the square of the matrix is computed. In the presence of complex 
elements the Hermitian transpose is performed. Finally, its reciprocal is calculated and 
the values obtained are translated to dB scale. The link overall estimated SNR in dB is 


then sent into the adaptive rate control block. 
9, Adaptive Rate Control 


In this block the instantaneous estimated channel SNR is compared with 
preloaded SNR values, which defines the limits of operation for each of the used seven 
modulation schemes. The circuit outputs an integer corresponding to the modulation 
scheme to be used by the transmitter and receiver units. The objective is to obtain an 


overall system with an error free data link. The output of this block is fed to the FEC and 
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Modulator bank and to the Demodulator and FEC bank to produce the adaptive rate 
control of the overall system. Also, its output is provided to the bit error rate block to 
correctly calculate the BER, since each modulation scheme has a different data block 


dimension. The values of the integer, designating rate ID, are presented in Table 9. 


Table 9. Rate ID Encodings, from [10]. 


Modulation scheme 




















eee an overall rate 
0 BPSK r=1/2 
1 QPSK r=1/2 
2 QPSK r=3/4 
3 16QAM r=1/2 
4 16QAM r=3/4 
5 64QAM r=2/3 
6 64QAM r=3/4 


C. MISO OFDM MODEL 


The MISO OFDM model is depicted in Figure 26. To avoid unnecessary 
description of blocks common to this model and the previous model, in the following 
section only the blocks that contain differences or are new in the MISO model will be 


described. 


4] 
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Figure 26. MISO OFDM Model. 
1. Space-Time Diversity Encoder 


The Space-Time Diversity Encoder is presented in Figure 27. This block receives 
the signal from the IFFT Input Packing block and produces two signals that are fed to the 
OFDM transmitters. The block itself is an embedded MATLAB® code S-function. It 
contains MATLAB® code that implements the Alamouti scheme described in Chapter II. 
The code is available in the Appendix. 
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Figure 27. Space-Time Diversity Encoder. 


2. OFDM Transmitters 


Both OFDM transmitters have the same configuration as the SISO OFDM 
transmitter. However, in the MISO model the transmitters are supplied with two different 
preambles, in contrast to the SISO model where the transmitter was supplied with a 
single short preamble. The first transmitter is supplied with the even preamble and the 


second transmitter is supplied with the odd preamble. 
3. MISO Fading Channel 


This block, shown in Figure 28, simulates a MISO multipath fading channel with 
two inputs and one output. As discussed in the previous chapter it is a particular case of a 
MIMO channel obtained by simply disconnecting one of the outputs. The values for the 
phases were selected randomly. The correlation parameters simulated by the gain 
elements were selected to be 0.5 to simulate an intermediate level of correlation between 


the antenna elements. Considering the used carrier frequency f, = 2.3GHz, the spacing 


d between antennas can be calculated using the relation [21] 


22d 
py 13( 74), (3.3) 


where J, (x) is the Bessel function of the first kind of order zero and / is the signal 


wavelength. Is this case the antennas spacing is d ~ 0.182. 
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The Rayleigh Fading elements simulate a multipath fading channel with 


parameters defined by the user. 
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Figure 28. MISO Fading Channel. 


4. Additive White Gaussian Noise Channel 


The AWGN channel is a replica of the AWGN component used in the Multipath 
Fading Channel with AWGN for the SISO model. Its sole function is to add noise to 
overall channel modeling as in the SISO model. The SNR in GB is defined by the user. 


5, Space-Time Diversity Combiner 


In this block, the received signal is detected after being adjusted by the channel 
estimation. In the first section, the channel is estimated by dividing the received preamble 
by the transmitted preamble. The channel’s estimations and the received signal with the 


DC component removed are then provided to the space-time combiner for signal 
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detection. This block is an embedded MATLAB® code S-function. It contains 
MATLAB® code that implements the Alamouti scheme at the receiver as described in 


chapter II. The code is available in the Appendix. 


Finally, the amplitude of the detected signal is adjusted by the Gain Compensator. 
The detected signal is divided by the sum of norm of the channels estimations to obtain 
the desired compensation. The detected signal is then sent to the Data Carries Extraction 


block for further processing similar to that of the SISO model. 
D. MIMO OFDM MODEL 


The MIMO OFDM model is shown in Figure 29. The blocks are similar to the 
previous cases apart from the MIMO Fading Channel and the Space-Time Diversity 


Combiner, which are presented next. 
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Figure 29. MIMO OFDM Model. 


1, MIMO Fading Channel 


This block, shown in Figure 30, simulates a MIMO multipath fading channel with 


two inputs and two outputs. The phase values were selected randomly. Again, correlation 


parameters simulated by the gain elements were selected to be 0.5 to simulate an 


intermediate level of correlation between the antenna elements. The Rayleigh Fading 


elements simulate a multipath fading channel where the delay and gain parameters are 


defined by the user. 
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Space-Time Diversity Combiner 
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Figure 30. MIMO Fading Channel. 


The Space-Time Diversity Combiner block, shown in Figure 31, as in the MISO 


model, detects the received signal as compensated by the channel estimation. First, the 


four possible channels are estimated by dividing the received preamble by the transmitted 


preamble. The channels’ estimations and the received signals with the DC component 


removed are then provided to the space-time combiner for signal detection. This block is 


an embedded MATLAB® code S-function. It contains MATLAB® code that implements 


the Alamouti scheme at the receiver as described in chapter II, for two receiving 


antennas. The code is available in the Appendix. 


Lastly, the detected signal is compensated in gain by the Gain Compensator 


depicted in Figure 32. The detected signal is divided by the sum of norm of the channels’ 
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estimations to obtain the desired compensation. The detected signal is then sent to the 
Data Carriers’ Extraction block for further processing analogous to that of the SISO and 


MISO models. 
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E. ITU CHANNEL MODELS 


For simulation purposes, the multipath channel models chosen for the executed 
simulations were those specified by the ITU. These channel models provide a variety of 
situations considered typical. Three user locations are considered: indoor, pedestrian and 
vehicular. The indoor location user is a fixed subscriber, thus its Doppler spread is null. 
The pedestrian user has normally a speed up to 3km/h, which will induce a maximum 
Doppler shift of ff, =6.389Hz, considering the use of a carrier frequency of 
f, =2.3GHz. The vehicular user is considered to have a speed between 60km/h and 
120km/h, which corresponds to Doppler shifts of f, =127.77Hz and f, =255.56Hz, 
respectively, for the same carrier frequency. For each one of these types of user, the ITU 
specified two profiles of multipath: A and B. Profile A has shorter delay spread when 
compared to profile B. Profile A replicates rural macro-cellular surroundings, while 


profile B reproduces an urban macro-cellular environment. For scenarios of micro-cells 
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with radius less than 500m profile A is also suggested. The values of delay and relative 
power for each of the users and profiles are presented in tables Table 10, Table 11 and 


Table 12 


Table 10. ITU Multipath Channel Models — Indoor, after [1]. 




















Tap Number _ Delay (ns) Saar Delay (ns) esa 
Channel A Channel B 

1 0 0 ; 5 

2 50 3 100 36 
3 110 -10 200 42 
= 170 -18 300 -10.8 
5 290 -26 500 -18.0 
6 310 32 700 ‘955 


Table 11. ITU Multipath Channel Models — Pedestrian, after [1]. 

















Tap Number __ Delay (ns) See ai) Delay (ns) Cass 
Channel A Channel B 

1 0 0 0 0 

2 110 -9.7 200 -0.9 
3 190 -19.2 800 -4.9 
4 410 -22.8 1200 -8.0 
5 2300 -7.8 
6 3700 -23.9 
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Table 12. ITU Multipath Channel Models — Vehicular, after [1]. 

















Tap Number _ Delay (ns) Rare Delay (ns) an) 
Channel A Channel B 

1 0 0 0 -2.5 
2 310 -1 300 0 

B) 710 -9 8900 -12.8 
4 1090 -10 12900 -10.0 
5 1730 -15 17100 -25.2 
6 2510 -20 20000 -16.0 


F. SUMMARY 


In the present chapter we have introduced the SISO, MISO and MIMO OFDM 
models used in this thesis to obtain the results presented later. A brief presentation on the 
ITU channel models used in the simulations was also provided. In the next chapter, we 


will address the results obtained from the several performed simulations. 
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IV. SIMULATIONS AND RESULTS 


A. INTRODUCTION 


In the previous chapter, a description of the several communication systems and 
multipath channel models used in this thesis were presented. In the present chapter the 
simulation settings used and the results obtained will be addressed. On the results, first 
the performance of the communication systems is discussed with respect to BER versus 


SNR and later the measured system’s capacity is presented. 
B. SIMULATION SETTINGS 


Table 13. Simulation Settings for Systems Benchmarking. 














Parameter Value 
Channel bandwidth B 3.5MHz 
Number of subcarriers L,,; 256 
Carrier frequency f. 2.3GHz 
Ratio of cyclic prefix time to useful symbol time G =T, / T, 1/8 
: ‘ Pr =9.5 
MISO & MIMO fading correlations 
Pp =9.5 





¢ =1.67; ¢, =2.11 


MISO & MIMO rand h 
S random phases @, = 0.25; g, = 0.89 





Number of data bits transmitted 10° 


The simulation settings were chosen to be those of the TEEE® 802.16-2004, also 
known as fixed WiMAX. Table 13 presents the significant parameter values of the 
OFDM systems and also the MISO and MIMO channel factors. Each system was tested 


under several ITU user-channel profiles presented in the previous chapter. The AWGN 
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SNR value was varied from OdB up to 44dB, in 1dB steps, for each modulation scheme 
presented in Table 9. For each SNR value 10° data bits were transmitted to measure BER 
values as low as 10°. The curves obtained were considered to be the system’s 
benchmarks under the several user-channel conditions. From these performance curves, 
after observing the values obtained for each modulation at a BER of 10“, we defined the 
SNR threshold vectors for the rate control. These values were used in the second part of 
the simulations to obtain performance figures of all systems with partial CSI feedback. 


The simulation settings in this phase are presented in Table 13. 
C; PERFORMANCE RESULTS 


In this section, the overall performance in terms of measured BER versus the link 


overall SNR is discussed for several user profiles and channel profiles. 
1; AWGN Channel Performance 


The systems were first tested with an AWGN channel model. These results 
provide a benchmark for comparison when the multipath channel effect is added to the 
simulations. The results obtained for this case are depicted in Figure 33 and Figure 34. 
The MISO system showed very little improvement when compared to the SISO system. 
On average, a gain less than 1/2dB was achieved. In the particular case of the BPSK 
signal, the MISO system actually performed worse than the SISO system. However, the 
MIMO system performed significantly better than both the SISO and MISO systems. On 


average a gain of 2dB was obtained. 
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Figure 33. Performance in AWGN Channel for PSK Signals. 
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Figure 34. Performance in AWGN Channel for QAM Signals. 
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2. AWGN plus Multipath Channel Performance 


In this general channel scenario, all ITU profiles presented in the previous chapter 


were simulated. In the next sections the relevant results are discussed. 
a. Indoor Channel A 


In this scenario, the results obtained were encouraging. With BPSK and 
QPSK the MISO system performed on average 1dB better than the SISO system, while 
the MIMO system showed a performance 3dB better than the MISO system. These 
results are presented in Figure 35. Figure 36 shows the results for the QAM signal 
families. In this case the MISO system shows a 2dB average performance improvement 
from the SISO system. Also, the MIMO system has an average performance 
improvement of 3dB when compared to the MISO system, and a 5dB improvement 


when compared to the SISO system. 
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Figure 35. Performance in AWGN plus Multipath Indoor A for PSK Signals. 
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Figure 36. Performance in AWGN & Multipath Indoor A for QAM Signals. 


b. Indoor Channel B 


In this simulation profile some significant results were obtained. Recall 
that the profile of channel B has a bigger time delay spread than the profile of channel A, 
more than twice to be more quantitative. This factor plays a big role in the systems’ 
performances. Observing both Figure 37 and Figure 38, the performance of the SISO 
system is not satisfactory, in particular for the case of QPSK r=3/4and QAM signals. 
Where comparison is possible, a 13dB average improvement was obtained when 
comparing MISO to SISO system. The MIMO system showed a_ performance 
improvement compared to the MISO system that varies between 3dB and S5dB, 
depending on the type of signal in use. Overall, in this particular user-channel profile, the 
MIMO system reached a peak performance improvement of 18dB when compared to the 


SISO system. 
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Figure 37. Performance in AWGN & Multipath Indoor B for PSK Signals. 
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Figure 38. Performance in AWGN & Multipath Indoor B for QAM Signals. 
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Cc. Pedestrian Channel A 


In the pedestrian profile, two different situations were considered: a 
moving and a stationary person. In this latter case, the MISO system performed worse 
than the SISO system for almost every type of signal. A 1/2dBto 1dB degradation was 
measured. On the other hand, the MIMO system performed better than the SISO system 


on average 2dB to 3dB for all types of signals. These results are depicted in Figure 39 
and Figure 40. 


For the case of a fast moving pedestrian, the SISO system performance 
was poor. Some BER curves showed unpredictable behavior, leading to its exclusion of 
use under this user-channel profile. Where it was possible to compare the MISO system 
performed better than the SISO system, except for the case of the 16QAM r=1/2 signal. 
Comparing the MIMO to the MISO system, the measured average improvement of 


performance was on the order of 5dB. 
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Figure 39. Performance in AWGN & Multipath Stopped Pedestrian A for PSK Signals. 
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Figure 40. Performance in AWGN & Multipath Stopped Pedestrian A for QAM Signals. 
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Figure 41. Performance in AWGN & Multipath Active Pedestrian A for PSK Signals. 
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Figure 42. Performance in AWGN & Multipath Active Pedestrian A for QAM Signals. 


d. Pedestrian Channel B 


Using the same methodology as in the previous section, simulations for 
both active and stationary pedestrians were carried out. For both cases the SISO system 
was unable to deal with multipath. Only the BPSK and QPSK r=1/2 signals were able 
to combat multipath. Because of this fact, no comparison of the SISO system to the other 
systems will be addressed. Also, in the case of the MISO and MIMO systems, the 
utilization of 64QAM signal appears to be impossible. As soon as the noise effect on the 
channel is negligible the multipath effect becomes clear, since the performance curves 
present horizontal asymptotes. Where comparison was possible, 2dB to 3dB average 
performance improvement was measured, comparing the MIMO to the MISO system in 


the stationary pedestrian case. In the active pedestrian situation 3dB to 4dB 


improvement was achieved. 
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Figure 43. Performance in AWGN & Multipath Stopped Pedestrian B for PSK Signals. 
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Figure 44. Performance in AWGN & Multipath Stopped Pedestrian B for QAM Signals. 
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Figure 45. Performance in AWGN & Multipath Active Pedestrian B for PSK Signals. 
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Figure 46. Performance in AWGN & Multipath Active Pedestrian B for QAM Signals. 


63 








e. Vehicular Channel A 


The IEEE® 802.16-2004 standard was not proposed for communication 
links under mobility, i.e., vehicular use. However, simulations under these conditions 
were performed to get a sense of the effects and to eventually reflect on solutions to 
combat the negative consequences. With this perspective in mind, we will only present 
two extreme cases under mobility, among all that were simulated. In this section the 


performance of the link under the user-channel vehicular A profile with 60km/h is 


addressed. In the following section the user-channel vehicular B profile with 120km/h is 


presented. 
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Figure 47. Performance in AWGN & Multipath Vehicular A for PSK Signals. 


The results obtained for the user-channel vehicular A profile with 
60km/h are somewhat promising in the MIMO case. The SISO and MISO systems 


clearly showed overall poor performance. In the cases where the systems could be 


compared, the MIMO system performed 3dB to 4dB better than the MISO system. For a 
64 


BER on the order of 10“, the MIMO system can be considered a viable solution, since 


only with the 16QAM r= 2/3 and 64QAM cases did it show poor performance. 
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Figure 48. Performance in AWGN & Multipath Vehicular A for QAM Signals. 


f Vehicular Channel B 


In this section, we present the results obtained under the user-channel 
vehicular B profile with speed 120km/h. It is clear that all systems performed poorly 
since none of them can combat the multipath and Doppler spread combined effect of this 
kind of channel. Also visible is the fact that the MIMO outperforms the MISO which 
outperforms the SISO system. Looking at the point where the noise effect on the channel 
is negligible (around 25dB SNR), the best achieved BER are 0.006, 0.02 and 0.04 for 
the MIMO, MISO and SISO systems, respectively, using the BPSK signal. This has 
tremendous effects on the channel throughput and confirms the fact that the IEEE® 


802.16-2004 standard has not been designed for mobile applications. 
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Figure 49. Performance in AWGN & Multipath Vehicular B for PSK Signals. 
Performance AWGN & Multipath Vehicular B 120km/h 
1.0E+00 

















BER 





























SNR (dB) 








—e— 16QAM 1/2 SISO 
—s— 16QAM 1/2 MISO 
—s— 16QAM 1/2 MIMO 
—x— 16QAM 3/4 SISO 
—*— 16QAM 3/4 MISO 
—s— 16QAM 3/4 MIMO 
—— 64QAM 2/3 SISO 
—— 64QAM 2/3 MISO 
— 64QAM 2/3 MIMO 
—e— 64QAM 3/4 SISO 
—s— 64QAM 3/4 MISO 
—a— 64QAM 3/4 MIMO 














Figure 50. 


Performance in AWGN & Multipath Vehicular B for QAM Signals. 
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3s AWGN plus Multipath Channel Performance with Partial CSI 
Feedback 


Table 14. SNR Threshold Vectors for Auto Rate Control Mode. 












































User-Channel Profiles System Vector 
SISO [12.5,17,18.5, 23.5,30.5,35.5] 
Indoor A MISO [12,16.5,17, 22.5, 26.5, 30.5] 
MIMO [8,12,13,18, 23, 27] 
SISO [ 25,-,28,-,-,-] 
Indoor B MISO [11.5,—,16, 24, 26.5, 33] 
MIMO [9,12°5,13,18;22,97] 
SISO [13.5,18,19.5, 25,33, 38] 
Pedestrian A not moving MISO [14,19, 20, 26, 31,36] 
MIMO [11,-,16,22,27,31] 
SISO [-,-,21,-,34, 44] 
Pedestrian A moving MISO [14,19, 20, 26,31, 36] 
MIMO [11,-,16, 22, 27,31] 
SISO [28, -,36,-,-,—] 
Pedestrian B not moving MISO [14,-,16,-,-,- 
MIMO [9,-,14, 24,-,-] 
SISO [28,-,36,-,-,-] 
Pedestrian B moving MISO [14,- 635 -| 
MIMO [9,-,14,24,-,—] 


In this phase of simulations, the systems were set in auto rate control mode. That 
is, using the performance curves obtained in the previous simulations, a set of SNR 
threshold vectors was defined and its values introduced under each system’s settings. 
These vectors, presented in Table 14, were defined for each of the ITU user-channel 
profiles and aim for a maximum BER of 10~. Taking into account the performances 


obtained in the previous simulations, the SNR threshold vectors were only defined for 
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indoor and pedestrian user-channel profiles. The simple AWGN channel was ignored 
since it is not realistic for indoor or pedestrian environments. Some vectors do not present 
values in some indices, meaning that a particular or several modulations (rate ID) can not 
be employed. If the missing value was located before an existing value, for simulation 
purposes it assumes the upper existing value. If there are one or more none existing 
values in the right side of the vector, for simulation purposes they assume extremely high 
values, so that the rate control block never selects a rate ID that would induce a BER 
above the desired value. Simulations for all the systems and user-channel profiles 
presented in Table 14 were executed to validate the defined vectors. In some few cases, 
an adjustment was made in some values to obtain the desired BER of 10~*. The vectors 
presented in Table 14 are the tuned values. Running simulations with the auto rate 
control mode enabled showed effective and as depicted in Figure 51, Figure 52 and 
Figure 53, the overall system’s performance was limited by the BPSK signal, as 
expected. In Figure 52 and Figure 53 it is also visible the degradation effect caused by the 
Doppler spread induced by an active pedestrian when compared to a stationary 


pedestrian. 
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Performance in AWGN & Multipath Indoor using CSI feedback 
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Figure 51. Performance in AWGN & Multipath Indoor using CSI Feedback. 
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Performance in AWGN & Multipath Pedestrian A using CSI feedback 
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Figure 52. Performance in AWGN & Multipath Pedestrian A using CSI Feedback. 
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Performance in AWGN & Multipath Pedestrian B using CSI feedback 
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Figure 53. Performance in AWGN & Multipath Pedestrian B using CSI Feedback. 


4. AWGN plus Multipath and Shadowing Channel Performance with 
Partial CSI Feedback 


One further step was taken to test the systems’ performance under adverse 
conditions. Simulations on all systems under AWGN plus multipath and shadowing were 
executed. The simulation settings were those presented in Table 13, plus a shadow 
standard deviation o, =8dB. The results for both indoor user-channel profiles are 
present in Figure 54. Again, the systems overall performance was limited by the use of 


BPSK modulation. Also visible is the auto rate control mode effect. Notice that errors 
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were measured beyond the BPSK curve, related to the other modulation schemes in use, 


but as expected below the established BER threshold of 10~*. 





Performance in AWGN & Multipath plus Shadowing Indoor using CSI 
feedback 
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Figure 54. Performance in AWGN & Multipath plus Shadowing Indoor. 


5. Achievable Data Rates 


The maximum data rates per modulation scheme can be calculated using the 
parameters presented in Table 13 and the parameter definitions available in [10]. The 


computed values are presented in Table 15. 
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Table 15. Maximum Data Rates per Modulation Scheme. 


. Maximum Channel Maximum Uncoded 
Modulation 


Scheme Data Rate Data Rate 

(Mbps) (Mbps) 
i lace 2.667 1.333 
paieeieaies 5333 2.667 
ee eat 5333 4.000 
dae 10.667 5.333 
dares 10.667 8.000 
i iets Aa 16.000 10.667 
Nt 16.000 12.000 


Table 16. Maximum Data Rates per User-Channel Profile. 


Maximum Channel Maximum Uncoded 





















































User-Channel Profiles System Data Rate Data Rate 
(Mbps) (Mbps) 
SISO 16.000 12.000 
Indoor A MISO 16.000 12.000 
MIMO 16.000 12.000 
SISO 10.067 335 
Indoor B MISO 16.000 12.000 
MIMO 16.000 12.000 
SISO 16.000 12.000 
Pedestrian A not moving MISO 16.000 12.000 
MIMO 16.000 12.000 
SISO 16.000 12.000 
Pedestrian A moving MISO 16.000 12.000 
MIMO 16.000 12.000 
SISO 10.067 D3395 
Pedestrian B not moving MISO 10.067 5.333 
MIMO 10.067 8.000 
SISO 10.067 5.550 
Pedestrian B moving MISO 10.067 5:335 
MIMO 10.067 8.000 
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From the values present in Table 15 and the taking into account the highest 
modulation scheme employed in each user-channel profile, the maximum achievable data 
rates per system and user-channel profile can be predicted. These values are presented in 
Table 16. Notice that these values are peak values and only consider a simplex link. If we 


consider a more realistic situation of systems using time division duplexing (TDD) with a 


down-link to up-link ratio of one to one (1) , the values are around half of those 


presented. Furthermore, the data throughput will be even less because of packet 
overhead. Nevertheless, it is visible from the values presented that the MIMO system 
capacity is superior to both the MISO and SISO systems capacities. In adverse situations, 
like larger delay spread multipath channels and/or presence of Doppler shift, the MIMO 
system can provide an increase of 50% or more in data rate when compared to the MISO 


and SISO systems. 
D. SUMMARY 


In the present chapter, we addressed the significant results obtained during this 
thesis work. In general, the MISO system outperformed the SISO system and was 
outperformed by the MIMO system under nearly all user-channel profiles. Furthermore, 
when the systems were set to use the auto rate control mode, based on the overall 
channel’s measured SNR, all adapted to the channel conditions. Cases of AWGN plus 
multipath and AWGN plus multipath with shadowing were simulated. In all, the systems 
reacted as expected, changing the rate ID as appropriate. Finally, predicted data rates for 
each system and user-channel profiles were presented as a measurement for the systems’ 
capacities. In the following chapter, we present a summary of the work performed during 


this thesis, the main results, conclusions, and suggestions for future work. 
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Vv. CONCLUSIONS 


The objective of the present thesis was to evaluate the MIMO OFDM 
performance and to reach its optimal data transmission. This was accomplished by 
selecting an OFDM standard and evaluating its performance under several user-channel 
profiles. The IEEE® 802.16-2004 standard was selected since it is a broad standard that 
permits several settings and implementation options, and also because it is an active area 
of research where its technical and commercial implementation is ongoing. The results 
presented in this thesis were obtained by systematic measurements using MathWorks'™ 


SIMULINK® R2008a software package models. 
A. SUMMARY OF THE WORK DONE 


A comprehensive study of MIMO OFDM material and wireless channel behavior 
was presented in Chapter II. The principles of OFDM were studied and its modern digital 
implementation was addressed. On the wireless channels concepts of path loss, fading 
and shadowing were introduced in order to understand the various effects that produce 
degradation in the transmitted signal. A broad introduction to MIMO systems was given 
starting from the basic SISO system, and passing through the intermediate MISO and 
SIMO systems. Following, the STBC technique known as Alamouti’s scheme was 


discussed. 


In total, three models were available by the end of the work. The first one was 
provided by the MathWorks™, and the remaining two were improved and modified 
versions. For simulation purposes, each system was complemented with MATLAB® code 
that enabled each simulation desired settings. In total, each SIMULINK® model has three 
corresponding MATLAB*® control programs. The simulations were performed using the 
ITU user-channel profiles, since they are the most frequently used power-delay profiles 
in simulations. The results obtained from the several simulations were presented in 


performance curves of BER vs. SNR, enabling comparison among the systems. 
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SIGNIFICANT RESULTS AND CONCLUSIONS 


Several significant results can be taken from this thesis: 


An increase in the number of antennas used in a wireless communications system 
enhances its performance and capacity. In this thesis, in most scenarios, the 2x1 
MISO system outperformed the SISO system and was outperformed by the 2x2 
MIMO system. 


User-channel characteristics under which wireless communications is tested or 
used have significant impact on the systems overall performance. In our case, 
using the IEEE® 802.16-2004 standard, it became clear that channels with larger 
delay spread are a bigger challenge to any system. The MIMO system proved its 
effectiveness in combating the multipath effect on the channels. The MISO 
system, under smaller delay spread conditions was also effective, but always less 
than the MIMO system. The SISO system experienced difficulties in combating 


the multipath effect in both large and small delay spread conditions. 


Also under the user-channel profile, the Doppler spread plays an important role. 
In indoor or pedestrian user profiles all systems were able to achieve a certain 
level of performance, where again the MIMO system was the best followed by the 
MISO system which was better than the SISO system. In the medium speed 
vehicular profile, only the MIMO and MISO systems were capable of 
counteracting the combined multipath and Doppler spread effect. Despite that, 
both systems were essentially limited to the use of the PSK family of signals. All 
systems, with the settings used, were incapable of dealing with high speed 


vehicular user-channel profiles. 


Measuring the overall SNR on a wireless communications link, as described in 
Chapter III, is a simple technique to obtain CSI. This enables only partial CSI and 
has limitations when auto rate control is desired. In this thesis this scheme proved 
effective while the noise present in the channel was considerable. When it became 
negligible and systems were left under the multipath effect alone, it was clear that 


the systems were restrained from reaching higher performances. 
76 


C. 


SUGGESTIONS FOR FUTURE WORK 


After the conclusion of this thesis, several areas for further research and work are 


indentified: 


Further develop the systems by increasing the number of antennas at either the 
transmitter or at the receiver, using the Alamouti’s scheme a 4x2 MIMO system 
with stacked STBCs. This system can be implemented using two 2x1 MISO 
systems with Alamouti’s scheme [1]. Eventually, a 4x4 MIMO system with 
stacked STBCs using two 2x2 MISO systems with Alamouti’s scheme may be 


realizable. 


Implement on the developed systems a more sophisticated CSI feedback that will 
enable the systems to achieve higher performances. The IEEE® 802.16e-2005 
standard, which extends the previous standard to allow usage under mobility, 
points to two techniques: codebook based feedback and quantized channel 
feedback [22]. In both cases the standard does not point out how to implement 
such feature, however two ways of doing so are: maximization of sum capacity 


and minimization of mean squared error [1]. 


In the developed system, we have only focused on the IEEE® 802.16-2004 
standard specifications. This led to the use of a specific bandwidth, number of 
subcarriers (FFT size), and carrier frequency throughout the work. Further work 
using the IEEE® 802.16e-2005 standard, where different bandwidths, carrier 
frequencies and number of subcarriers are allowed, might show performances 


under certain conditions that are not yet understandable. 


The EC department presently has IEEE® 802.16-2004 standard equipment with 
one base station antenna and two subscribers, each with a single antenna. It is 
suggested to further invest in this equipment and the execution of field 
measurements to validate software models, but also to gather information to better 


tune future systems. 
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On a military perspective, it is becoming clearer every day that we will further 
and further make use of broadband wireless communication systems. In the naval 
environment, the VHF and UHF bands which naval forces use for tactical 
communications are stressed and overloaded. Interference is a constant factor 
when operating in a task force. The solution will be broadband wireless 
communication systems. A good example is the American Digital Wideband 
Transmission System (DWTS). This system has high capacity, can be used with 
or without security features, is LOS, ship-to- ship or ship-to-shore, and operates in 
the upper UHF band (1350-1850MHz), with data rates up to 2.3 Mbps [23]. Most 
of the time naval task forces sail under formations or screens to enable 
optimization of the combined weapon systems. Under such conditions multipath 
effects, such as those we have studied in this thesis are a reality. The results 
presented in this work are a good insight to what can happen to tactical naval 
communications using the DWTS or similar systems. This author has no 
knowledge of multipath channel models for naval environments. It would be 
desirable to collect wireless channel measurements that permit the design of 
models applicable to naval task forces, and from there study its effects on 


broadband communications to further improve existing and future systems. 
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APPENDIX S-FUNCTIONS CODE 


MISO 





Encoder 


function [fantl, ant2] = stbcenc(u) 
% STBCENC Space-Time Block Encoder 
% Outputs the Space-Time block encoded signal per antenna. 














complex (zeros (size(u))); 
ant2 = antl; 





% Alamouti Space-Time Block Encoder, G2, full rate 
% G2 = [sO sl; -sl* s0*] 
for i = l:size(u,2)/2 
sO = u(:, 2%*1i-1); sl = u(:, 2%*1); 
antl(:, [2*1i-1 2*1i]) [sO -conj(s1)]; 
ant2(:, [2*i-1 2%*i]) [sl conj(s0)]; 
end 


Decoder 








function z = stbcdec(chEstl, rx, chEst2) 
% STBCDEC Space-Time Block Combiner 

















N 2; M= 1; 
Zz = complex(zeros(size(rx))); 
zO = compl ees renee A M)); zl = 20; 
% Space Time Combiner 
for i = l:size(rx,2)/2 
zO(:, M) = rx(:, 2*i-1).* conj(chEst1l(:, 2*i-1)) + 
con] (t#x(s, 2*2)).* chEst2 (3, 2*1);7 
zZzl(:, M) = rx(:, 2*i-1).* conj(chEst2(:, 2*1i-1)) - 
cony(rx(:, 2*2)).* chEstl(:, 2*21); 
z(s, [2*i-1 24*2]) = [20 z1] 
end 
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MIMO 
























































chi 





ESt4) 


Encoder 
function [fantl, ant2] = stbcenc(u) 
% STBCENC Space-Time Block Encoder 
% Outputs the Space-Time block encoded signal per antenna. 
N = 2; 
antl = complex(zeros(size(u))); 
ant2 = antl; 
% Alamouti Space-Time Block Encoder, G4, full rate 
% G4 = [s0 sl; fas s0*;s0 sl;-—sl1* s0*] 
for i = l:size( )/2 
sO = u(:, oe Ae sl = u(:, 2*1); 
antl(:, [2*i-1 2*i]) = [sO -conj(s1)]; 
ant2(:, [2*i-1 2%*i]) = [sl conj(s0)]; 
end 
Decoder 
function z = stbcdec(chEstl, chEst2, rxl, rx2, chEst3, 
% STBCDEC Space-Time Block Combiner 
N 2; M= 1; 
Zz = complex(zeros(size(rxl))); 
zO = complex(zeros(size(rxl,1), M)); zl = 20; 
% Space Time Combiner 
for i = l:size(rxl1,2)/2 
z0 ( M) = rxl(:, 2*1i-1).* conj(chEst1l(:, 2*i-1)) 
econ (reLcs, 2°2)).* cHEStA (ey, 2*i) F 24 
LE2Cl,. 2*1-1).% con) (chkst3:(s,. 2*i-1)) 
conj(rx2(:, 2*i)).* chEst4(:, 2*i); 
Zl(:, M) = rxl(:, 2*i-1).* conj(chEst2(:, 2*1i-1)) 
con] (2x1 Cr, 2%*2)).* GhESti-(s,. 242) + ws 
rx2(:, 2*i-1).* conj(chEst4(:, 2*i-1)) 
conj(rx2(:, 2*i)).* chEst3(:, 2*i); 
z(:, [2*i-1 2*i]) = [z0 zl]; 
end 
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